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UNIFIED  THEORY  OF  THE  POWER  SPECTRUM  OF  LONG  WAVELENGTH 
IONOSPHERIC  ELECTRON  DENSITY  IRREGULARITIES 

INTRODUCTION 


The  study  of  electron  density  irregularities  and  structures  in  the 

ionosphere  has  divided  itself  naturally  into  three  regimes:  the  low,  middle 

and  high  geomagnetic  latitudes.  These  regimes  have  fundamentally  different 

sources  for  the  irregularities,  associated  with  the  degree  of  coupling  to 

higher  altitude  magnetospheric  phenomena.  The  coupling  is  related  to  the 

orientation  of  the  geomagnetic  fields  as  a  function  of  magnetic  latitude.  At 

low  latitudes,  where  the  geomagnetic  field  tends  to  be  horizontal,  coupling  to 

higher  altitude  magnetospheric  disturbances  is  inhibited.  At  high  latitudes 

the  more  vertical  magnetic  field  promotes  strong  magnetospheric-ionospheric 

coupling.  (For  recent  reviews  of  ionospheric  structures  and  irregularities 

see  Fejer  and  Kelley,  1980,  Ossakow ,198 1 ,  and  Keskinen  and  Ossakow,  1983.)  A 

fundamental  quantity,  which  can  be  measured  experimentally  and  computed 

theoretically,  characterizing  ionospheric  structures  and  turbulence  is  the 

power  spectrum  of  the  density  fluctuations.  A  theory  which  extends  the 

Kolmogorov  picture  of  fluid  turbulence  to  unstable  plasma  systems  [Sudan  and 

Pf Irsch, 1982]  has  recently  been  applied  to  Type  II  electron  density 

irregularities  in  the  equatorial  electro jet  [ Sudan  1983].  The  concepts 

underlying  this  treatment  are  a  development  of  previous  studies 

[Kulsrud  and  Sudan,  1981;  Sudan  and  Keskinen,  1977,  1979  Keskinen,  1981]. 

This  theory  [Sudan,  1983]  predicts  the  absolute  magnitude  of  the  power 

spectrum  as  a  function  of  wavelength  in  terms  of  a  strength  parameter 
2  2 

v^YQ/k  Cg  =  S  which  can  also  be  interpreted  as  a  Reynolds  number  that  defines 

the  level  of  the  turbulence;  v,  is  the  ion-neutral  collision  frequency,  C_  is 

1  S 

the  sound  speed,  k  »  2nr / A  .  X  is  the  largest  scale  size  involved  in  the 

o  o  o 

turbulence  and  Y0(^0)  Is  the  linear  growth  rate  of  the  fluctuation  whose 
wavelength  is  XQ  •  The  purpose  of  thi3  letter  is  to  show  that  this  theory 
applies  equally  well  and  almost  without  significant  modification  to  equatorial 
spread  F,  natural  auroral  electron  density  irregularities  in  the  E  and  F 
regions,  and  artificial  irregularities  in  barium  clouds  for  wavelengths  above 
the  ion  gyro-radius  and  below  say  -  10km  .  The  predictions  of  this  theory  are 
in  good  agreement  with  the  accumulated  radar  observations  [Farley,  1979; 
Hanuise  and  Crochet ,  1981a,  1981b]  and  the  more  recent  rocket  observations 
[Fejer  and  Kelley,  1980;  Kelley,  et.  al. ,  1982;  Prakash  et .  al. ,  1972  Keskinen 
et  al. ,  1981;  Singh  and  Szuszczewicz ,  1983;  Rodriguez  and  Szuszczewicz ,  1983]. 
Manuscript  approved  October  4,  1983. 
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THEORY 

For  low-frequency  nodes  che  .onospheric  plasma  in  the  F  and  F  region  can 

be  node  lied  by  two-fluid  equations  in  which  quasi-neutrality  is  assumed 

a  ■  a.  ■  a  ,  the  pressure  variations  are  taken  to  be  isothermal,  and  the 
e  1  o 

electron  and  ion  inertia  are  both  neglected.  Because  8nn(T  +  T.)  <<  B  “ 

e  1  o 

where  B^  is  che  Earth's  magnetic  field,  the  perturbed  magnetic  fluctuations 
vanish  and  the  electric  field  fluctuation  are  electrostatic  E  *  -7p  and  we  can 
write 
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In  this  set,  the  e  and  i  subscripts  refer  to  electrons  and  ions, 

respectively,  u,  is  the  charged  particle  neutral  collision  frequency,  g  is 
i ,  e  '*• 

the  acceleration  due  to  gravity  and  the  other  symbols  have  and  their 
conventional  meaning.  In  the  continuity  equation  (1)  V,  is  the  particle  drift 

across  the  magnetic  field.  In  the  E-region  where  ions  are  unmagnetized, 

_2 

v./n,  -  10,  and  electrons  are  magnetized,  v  1 2  -  3  x  10  2  (  a.  Is  the 

ii  e  e  i 

cyclotron  frequency),  the  electrons  drift  at  a  velocity  given  by 

2 

V.  •  cE  x  B  /B  where  E  is  the  ambient  vertical  static  field.  The  ions  are 
—a  •»o  ~o  o  — o 

stationary  and  the  electrons  drift,  V,  results  in  the  electrojet  current.  In 

®  » 3 

the  F-region,  on  the  other  hand,  v.  /£,  -  10  ,  v  /fl  -  10  i.e.  ,  both 

li  e  e 

electrons  and  ions  are  magnetized  and  the  ambient  field  Eq  causes  the 
electrons  and  ions  to  drift  at  the  same  rate,  so  that  no  electric  durrent  is 
generated  by  E  .  However  the  much  weaker  gravitational  field  causes  the  ions 

''O 

to  drift  with  respect  to  the  electrons.  Thus  in  the  electron  drift  frame  we 
miC  2 

set  ■  (— — )  £  x  Bq/Bq  the  ion  drift  velocity  and  V  -  in  Eq.  (I). 
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*  V 


Consider  an  electrostatic  wave  54  *  j  exp  i(k.x  -  ^t)  propagating  in 
this  medium  such  that  kL  »  where  L  is  the  scale  length  of  the  ambient 
electron  density.  Such  a  wave  will  be  Landau  damped  by  electrons  free- 


streaming  along  SQ  unless 


J  >  *  ,ve>  k  j 


)/3  and  v  *  (2T  /a  )  1/2  is 


o"  0  e  x  e  e' 

the  electron  thermal  velocity.  Mow  the  low-frequency  modes  under  discussion 

Thus 


have  Re^,  3  k*V ,  =  k  V,  . 
k  ~d  id 


k  / k, 

r  1 


vv. 
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,10  ,  E  region 

’ 10~^ ,  F  region 


It  is  clear  that  these  irregularities  have  an  almost  two-dimensional 
structure  in  both  the  E  and  F  regions  and  we  may  consider  the  turbulent 
interactions  between  the  different  modes  essentially  as  a  two-dimensional 
system  and  isotropic  in  the  plane  perpendicular  to  3 

At  the  equator  the  magnetic  field  has  no  dip  and  the  ambient  electron 

density  gradient  7nQ,  and  the  mean  electrojet  drift  are  orthogonal  to 

each  other.  This  special  geometry  has  been  utilized  in  the  study  of 

electro jet  driven  irregularities.  However  in  the  auroral  regions,  while 

V  ,  and  3  are  orthogonal,  7n  and  3  are  not  because  of  the  large  dip  of  the 
-a  ~o  0-0 

earth's  magnetic  field.  In  this  situation  because  the  modes  of  interest  still 
have  k^  <<  k^  ,  the  component  of  the  density  gradient  along  the  magnetic 
field  BQ.Vno/Bo  ,  plays  no  role  in  the  dynamical  evolution  of  these  modes  so 

long  as  (k«B  )(B  «7n  /n  B2)  >>  1  (we  do  not  consider  j  •  B  *  0  driven  modes 

— O  — O  OOO  ^  v 

here).  The  auroral  case  is  equivalent  to  the  equatorial  situation  if  7nQ  is 
replaced  by  V^no  which  is  the  component  of  the  gradient  perpendicular 

to  B  .  Thus  the  analysis  for  the  equatorial  S  and  F  regions  could  be 
extrapolated  to  any  latitude  provided  the  appropriate  density  gradient  is 
taken  into  consideration. 

In  the  set  (1)  -  (4)  the  principal  nonlinear  interaction  takes  place 
through  the  term  7*nV  in  Eq.  (1).  Denoting 

*k,m  ”  ^(I772  *XP  iUt  ‘  ^ 

\tUi  '  $75^  n(x~’  0  exp  1  (,J,t  ‘  ^ 


It  can  be  shown  that  this  system  of  aquations  can  be  written  in  the  fora 
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with  k>>  »  k  -  k'  and  m  j  -  S  •  In  the  equatorial  I-region  we  have  shown 
previously  [Sudan  and  Keskinen,  1977;  Sudan,  1983]  that  the  Hall  current 
driven  E  x  3  instability  has  frequency  and  growth  rate 


that 
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(6a) 

(6b) 


(6c) 


To  extend  these  results  to  the  E-region  at  arbitrary  latitude  one  need  only 

replace  3n  /3z  with  7  n  .  The  nonlinear  terms  remain  unchanged, 
o  i  o 

The  Rayleigh-Taylor  and  E  x  B  gradient  drift  instabilities  have  been 
invoked  to  explain  low  and  high  latitude  F-region  density  irregularities, 
respectively  (see  recent  reviews  of  Ossakow,  1931  for  the  low  latitude 
ionosphere  and  Keskinen  and  Ossakow,  1983  for  the  high  latitude  case).  We  now 
show  that  the  structure  of  eq.  (5)  is  unchanged  for  Rayleigh-Taylor 
and  E  x  B  gradient  drift  modes.  In  the  drift  frame  of  the  electrons  the  ion 
drift  is  given  by  [x,y,z  are  the  vertical  east/west,  and  north/south 
coordinates;  Bq  »  Bqz  at  the  magnetic  equator] 

£.  x  £/%  l°w  latl-tude 

^d  =  ~io  "  ^(v. /a. )(cE  /B  )  high  latitude  ^ 

11  o 

and  for  low  latitudes  the  perturbed  electron  and  ion  drifts  are 


.■:1 
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5s  ■  -i7.  B  k*V.  fn 
1  o  ■“  ~io _ 

•j.  c  U~  n 
1  o 

to  linear  order  in  5n  .  Furthermore  from  (9)  and  (3) 


7.5V.  -  -vj  (m^.)*1  7.(e75s  +  T.T5n) , 

"1  0 

and  (9)  (10)  and  (11)  when  substituted  into  (1)  yield  (5)  with 
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Note  that  the  expressions  for  the  nonlinear  matrix  element  *  of  (6c) 

and  (12c)  are  identical  in  structure.  Similar  relations  can~be  obtained  for 
the  high  laltitude  case. Thus  the  nature  of  the  nonlinear  interaction  of  the 
Rayleigh-Taylor  and  _E_  x  _3.  modes  in  the  equatorial  and  high  latitude  F-region 
is  identical  to  the  E_  x  _B_  instabilities  of  the  E-region. 

Equation  (5)  defines  the  interaction  of  modes  k,  k',  k"  .  If  Reu^  is 
non-dispersive  i.e.  3uj, /3k  is  independent  of  k  ,  all  possible  k'  and  k"  will 

K  ^  v  s' 

interact  with  k  provided  only  that  they  satisfy  the  triangular 

equality  k"  ■  k  -  k-  .  If  ^  is  dispersive  than  only  those  modes  which  also 
simultaneously  satisfy  -  u^.,  contribute  to  the  interaction. 

Depending  upon  the  shape^of  w(k)i  this  reduces  the  number  of  interacting  modes 
from  a  continuous  to  a  small  discrete  set.  The  physical  repercussion  of  this 

5 


difference  in  Che  number  of  interacting  modes  is  to  cause  a  strong  interaction 
in  the  non-aispersive  case  which  results  in  the  nonlinear  self-damping  of 
the  mode [Kadomtsev,  !9o3;  Kraichnan,  1959 1 .  This  self-damping  may  also  be 
interpreted  as  a  spread  ii*.  i  *  7.  in  the  eigenf requency  ^  .  For  the 

dispersive  case  the  interaction  is  weak  and  the  linear  eigenf requencies 
are  retained  even  in  the  nonlinear  development  (so  called  weak-turbulence,  see 
Sagdeev  and  Galeev,  1969).  For  the  non-dispersive  modes  under  consideration 
here  the  turbulence  is  considered  "strong". 

Because  we  are  dealing  with  fully  developed  turbulence,  fluctuations  with 

all  possible  wavelengths  in  a  certain  range  are  present.  Define  the  power 

spectrum  I,  ,  ■  <|n,  |2/n  2>  where  <  >  denote  an  ensemble  average, 

a)  k ,  u>  <o 

Also  1,^  *  "d'u  1^,  and  the  energy  in  the  fluctuations  'Cy,  =  where  u  is  a 
constant.  The  total  energy  density 

U  (x)  *  u_'d2k  I  k 

00 

*  2mu^’d(lnk)  k2!^  for  isotropic  power  spectrum 


Thus  s(k)  ■  k2I^  is  the  energy  in  the  logarithmic  interval  d  (Ink)  *  dk/k.  An 
estimate  for  Che  self-damping  from  a  dimensional  analysis  of  (5)  proceeds 
as  follows.  Assuming  thac  the  interaction  is  dominated  by  modes 
with  k  -  k'  -  Is."  , 


‘  k 


v 


-  V. 


k,  k' ,k- 


In^lkik  = 


3kVd(«<2 


(13) 


2  1/2 

where  we  have  employed  (6c)  set  4k  v  (  and  <|nkl2/n^>  -  kl£  .  Furthermore 
from  (2)  and  (4)  we  find  in  the  linear  approximation  that 


\  M  S(W  ^  x  2  nk 
“kT-~ 
o 

so  that 


(14) 


rk  5 


A(\ 
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where  AV  is  the  velocity  of  an  eddy  of  scale  size  x  .  The  self-damping  is 
therefore  related  to  the  inverse  of  the  "eddy-turnover"  time  r,  • 


Xolmogorov  5  theory  ot  fluid  turbulence  states  that  the  rluid  energy  contained 
in  the  interval  .ik  is  transferred  into  the  interval  in  an  eddy  turnover 


tine  t,  •  This  is  to  sav  chat  I'd,  dink.);,  7,  is  conserved  in  the  inertial 
regime.  In  plasma  turbulence  however  the  energy  nay  be  augmented  or 
diminished  by  the  linear  growth  or  damping  of  the  waves.  Thus  we  have 
[Kulsrud  and  Sudan,  1982;  Sudan  and  Ffirsch,  1982] 


d  C^V,)  *  v,  k-I  (16) 

dink  *  *  *  * 


where  y,^  is  Che  linear  growth/ damping  rate.  Substituting  for  7,  from  (13)  we 
obtain  a  first  order  equation  which  is  easily  solved.  Now  y,  =  y  k^G  is  valid 

rC  O 

for  both  E  and  F  regions  with  appropriate  choice  of  vq  and  G  [See  (6b) 
and(2b)].  Define  a  strength  parameter  S*y  /k2G  where  k0  is  the  longest 
wavelength  in  the  system.  The  solution  of  (5)  for  I.  *  0  is  then  given  by 


ICO  -  x*8/3:i-x-:/3  -  1  (:c4/3  -  l)-.2 


I(x)  -  3-k-  (koVd/Yo)2l(k) 


where  x  ■  k/k  •  Thus  the  spectrum  of  (17)  originally  established  for  the 
equatorial  S-region  [Sudan,  1983]  should  be  universally  valid  for  the  E  and  ? 
regions  at  all  latitudes  for  almost  two-dimensional  modes  governed  by  the  set 
of  Equations  (1)  to  (4).  Notice  that  the  difference  in  the  nature  of 
instabilities  does  not  reveal  itself  in  the  spectrum  because  the  nonlinear 

term  7»nV_  which  established  the  interaction  between  the  modes  is  the  same. 

v  V2 

For  x  >  S  the  modes  are  linearly  damped  but  the  spectrum  will  extend 

3/*» 

to  x  -  S  before  it  vanishes.  Then  the  range  of  damped  modes  excited  by 

1/-  1/  2 

energy  transfer  from  unstable  modes  is  given  by  ^x  -  S  ‘  -  S 

DISCUSSION 

For  parameters  typical  of  the  equatorial  and  high  latitude  E  and  F 
4  6 

regions  S  >10  -  10  and  \  ■  2y/k  may  range  from  1-10  km.  In  the  inertial 

®  -8/3  0 

range  therefore  I(k)  «  k  .  This  power  law  dependence  and  associated 
spectral  index  is  consistent  with  previous  theoretical  and  experimental 
studies.  Fig.  1  shows  a  sample  power  spectrum  of  the  density  fluctuations 
resulting  from  the  numerically  obtained  evolution  of  the  Rayleigh-Taylor 


instability  in  the  equatorial  7-ragion  ionosphere  [Keskinen  et  al. ,  1931;. 
Good  agreement  is  seen  between  the  theoretical  prescription  for  the  power 
spectrum,  eq.  (17),  and  the  numerical  results.  Similar  power  laws  and 
spectral  incites  are  also  observec  in  theoretical  and  numerical  studies  : :  the 
Z_  x  _3  in  the  auroral  ionosphere  [Keskinen  and  hssakow,  1982,  1983]  and  in 

plasma  (bariua)  clouds  [Keskinen  et  al].,  193CJ.  In  addition,  the 
theoretically  derived  power  spectrua,  eq.  (17),  is  consistent  with  rocket 
observations  of  interaediate  scale  size  density  fluctuations  in  equatorial 
spread  F  [Keskinen  et  al,  1981,  Szaszczewicz  et  al. ,  1980;  Kelley  at  al. , 
1982;  Kelley  et  al.,  1982].  Fig.  2  shows  a  saaple  experimentally  obtained 
power  spectrum  taken  from  Keskinen  et  al.  [1981].  Similar  power  laws  are  also 
observed  with  satellite  measurements  in  the  high  latitude  ?  region  ionosphere 
[F.  Szuszczewicz ,  private  communication] . 
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Numerically  obtained  vertical  vave  number  spectra  [Keskinen  et  al, 

1981]  of  Rayleigh- Taylor  instability.  The  symbol  k  *  2o/960m  and 

o 

the  cutoff  valve  k  *  30.  The  continuous  line  is  obtained  from  eq. 
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FAIRBANKS ,  AX  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  3ROWN  (UNCLASS  ONLY) 

gte  syl vania,  inc. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  C03LE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61320 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  3EAU REGARD  STREET 
ALEXANDRIA,  VA  22311 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  SAUER 
01CY  ATTN  HANS  UOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  4  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LI3RARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01CY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  3ARBARA,  CA  93102 
01CY  ATTN  DAS I AC 
0  ICY  ATTN  WARREN  S.  KNAPP 
0  ICY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
0  ICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  4  SPACE  CO.,  INC 
P.O.  3CX  504 
SUNNYVALE ,  CA  94088 
01CY  ATTN  DEPT  60-12 
01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  4  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01CY  ATTN  MARTIN  WALT  DEPT  52-12 
01CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01CY  ATTN  RICHARD  G.  JOHNSCN  DEPT  52-12 
01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 


Old 

ATTN 

S.  HARRIS 

Old 

ATTN 

J.  MOULE 

Old 

ATTN 

GEORGE  MROZ 

01CY 

ATTN 

W.  OLSON 

Old 

ATTN 

R.W.  HALPRIN 

Old 

ATTN 

TECHNICAL  LI3P.ARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CT  ATTN  P.  FISCHER 

ATTN  V.F.  CREVIER 
ATTN  STEVEN  L.  G'JTS CHE 
ATTN  R.  30GUSCH 
ATTN  R.  HENDRICK 
ATTN  RALPH  KILB 
ATTN  DAVE  SOULE 
ATTN  F.  FAJEN 
ATTN  M.  5CHEIBE 
ATTN  CONRAD  L.  LONGMIRE 
ATTN  3.  WHITE 


01CT 

01CY 

01CT 

01CT 

Old 

Old 

Old 

Old 

Old 

Old 


MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  SEW  MEXICO  37 106 
Old  R.  STELLINGWERF 
Old  M.  ALME 
Old  L.  WRIGHT 


MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

Old  ATTN  JOHN  MORGANSTSRN 
Old  ATTN  G.  HARDING 
Old  ATTN  C.E.  CALLAHAN 


MITRE  CORP 

WESTGATE  RESEARCH  PARX 
1820  DOLLT  MADISON  3LVD 
MCLEAN,  VA  22101 
Old  ATTN  W.  HALL 
Old  ATTN  W.  FOSTER 


PACIFIC-SIERRA  RESEARCH  COR? 
12340  SANTA  MONICA  3LVD. 

LOS  ANGELES,  CA  90025 

Old  ATTN  E.C.  FIELD,  JR. 
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PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
313  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LA3 


PHOTCMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01301 

0  ICY  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  WA  98009 
Old  ATTN  E.J.  FREMOUW 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 


R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
01CY  ATTN  FORREST  GILMORE 

ATTN  WILLIAM  B.  WRIGHT,  JR. 
ATTN  ROBERT  F.  LELEVIER 
ATTN  WILLIAM  J.  KARZAS 
ATTN  H.  ORY 
ATTN  C.  MACDONALD 
TURCO 
DeRAND 
TSAI 


01CY 

01CY 

01CY 

01C? 

Old 

01CY 

01CY 

Old 


ATTN 
ATTN  L. 
ATTN  W. 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 


RAYTHEON  CO. 

528  30ST0N  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 


RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01CY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  3URNS 
01CY  ATTN  C.  SMITH 
OICY  ATTN  R.  TSUNODA 
01CY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  J.  VICKREY 
OICY  ATTN  RAY  L.  LEADA3RAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  R.  LIVINGSTOS 
OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 

OICY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN’  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D/  STQCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
3URLINGT0N,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 

PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 
NOTED) 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

Dr.  ?.  MANGE  -  CODE  4101 
Dr.  P.  GOODMAN  -  CODE  4130 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  W.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  D.  FARLEY 
DR.  H.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  K.  MAEDA 
DR.  S.  CURTIS 
DR.  M.  DUB IN 

DR.  N.  MAYNARD  -  CODE  696 
COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 
ENVIRONMENTAL  LABORATORY 
30ULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  C.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  X.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 


INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW1 6-262 
CAMBRIDGE,  MA  02139 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  E.  FRIEMAN 


23 


STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  ?.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  ,ffl  ..  .  _ 

DR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AX  99701 
DR.  L.E.  LtE 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  37 54 5 
DR.  M.  PONGRAT2 
DR.  D.  SIMONS 
DR.  G.  3AXASCH 
DR.  L.  DUNCAN 
DR.  P.  BERNHARDT 
DR.  S.P.  GARY 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  S.  OTT 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
DR.  R.  GREENWALD 
DR.  C.  MENC 

UNIVERSITY  OF  PITTS3UP.GR 
PITTSBURGH,  PA  15213 
DR.  N.  Z/.3USKY 
DR.  M.  3I0NDI 
DR.  E.  OVERMAN 


UNIVERSITY  OF  TEXAS 
AT  DALLAS 

CENTER  FOR  RESEARCH  SCIENCES 
P.O.  30X  638 
RICHARDSON.  TX  7 5080 
DR.  R.  HEELI5 
DR.  y.  HANSON 

j.?..  McClure  % 

UTAH  STATE  UNIVERSITY 
4TH  AND  3TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  X.  BAKER 
DR.  R.  SCHUNK 
DR.  J.  ST. -MAURICE 

PHYSICAL  RESEARCH  LABORATORY 
PLASMA  PHYSICS  PROGRAMME 
AHMEDA3AD  380  0C9 
INDIA 

P.J.  PATHAK,  LIBRARIAN 

LABORATORY  FOR  PLASMA  AND 
FUSION  ENERCY  STUDIES 
UNIVERSITY  OF  MARYLAND 
COLLECE  PARK,  MD  20742 
JHAN  VARYA N  HILLMAN , 
REFERENCE  LIBRARIAN 
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